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ABSTRACT 

We present the analysis of the integrated spectral energy distribution (SED) from the ultraviolet (UV) to the far- 
infrared and Ha of a sample of 29 local systems and individual galaxies with infrared (IR) luminosities between 10 11 L© 
and 10 118 L© . We have combined new narrow-band Ha+[N n] and broad-band g, r optical imaging taken with the Nordic 
Optical Telescope (NOT), with archival GALEX , 2MASS, Spitzer , and Herschel data. The SEDs (photometry and 
integrated Ha flux) have been fitted with a modified version of the MAGPHYS code using stellar population synthesis 
models for the UV-near-IR range and thermal emission models for the IR emission taking into account the energy 
balance between the absorbed and re-emitted radiation. From the SED fits we derive the star-format ion histories 
(SFH) of these galaxies. For nearly half of them the star-formation rate appears to be approximately constant during 
the last few Gyrs. In the other half, the current star-format ion rate seems to be enhanced by a factor of 3-20 with 
respect to that occured ~1 Gyr ago. Objects with constant SFH tend to be more massive than starbursts and they are 
compatible with the expected properties of a main-sequence (M-S) galaxy. Likewise, the derived SFHs show that all 
our objects were M-S galaxies ~1 Gyr ago with stellar masses between 10 10 ' 1 and 10 11-5 M 0 . We also derived from 
our fits the average extinction (A v = 0.6 — 3 mag) and the polycyclic aromatic hydrocarbons (PAH) luminosity to Lir 
ratio (0.03-0.16). We combined the A v with the total IR and Ha luminosities into a diagram which can be used to 
identify objects with rapidly changing (increasing or decreasing) SFR during the last 100 Myr. 
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1. Introduction 

Galaxies with high infrared (IR) luminosities (Lir = L(8 — 
1000 pm) > 1 0 11 L 0 ) are rare in the local Universe (e.g., 
Le Floc’h et al. |2005 ), and yet they are a cosmologically 
important class of objects because they dominate the star- 
formation rate (SFR) density at high- 2 . They can be classi¬ 
fied as luminous (Lir = 10 11 —10 12 L 0 ) and ultra-luminous 
(Lir > 10 12 L e ) IR galaxies (LIRGs and ULIRGs, respec¬ 
tively). LIRGs domina te the SFR density at 2 
ULIRGs do so at 2 ^ 2 (Perez-Gonzalez et al. 


1, while 


et al. 


2007) 


2005 Caputi 


The bulk of the IR luminosity of U/LI RGs is pro¬ 


duced by st rong star-formation (SF) bursts (|Sanders & 
Mirabel , |1996), although, the number of U/LIRGs with an 
active galactic nuclei (AGN) detected in their optical spec¬ 
tra increases with the IR lumin osity , reaching ~50% when 


Lir > 10 12,3 L o (Yuan et al. 2010). Similarly, the relative 


^20-25% of the luminosit y of local ULIRGs (Farrah et al. 


2007; Nardini et al. 2009). 


These episodes of intense SF and high IR luminosi¬ 
ties are mainly triggered by major mergers invol ving gas- 
rich progenitors in the case of local ULIRGs (Sanders 


& Mirabel 1996). Actua lly a majority (>80 %) of local 


ULIRGs are mergers (e.g., Veilleux et al.|2002| ). This is not 


true for local LIRGs which have more varied morph ologies 


et al. 


- isolated disks, disturbed spirals, or mergers (e.g.,[Arr ibas 

rl 


2004) Alonso-Herrero et al.||2006 Hung 
Hammer et al. ~ p005p suggest that episodic SF bursts due 
to minor mergers^or gas infall) could enhance the IR lu¬ 
minosity above the LIRG threshold. Moreover, according 
to TV-body simulations, the most likely remna nts of these 


contribution of AGN to the bolometric luminosity increases 
with increasing IR luminosity provi ding less than 2 -15% of 
the total luminosity of local LIRGs ([Pereira-Santaella et al.l 

2011 Petrie et al., 2011 Alonso-Herrero et al.[ 2012a| ), and 


minor interactions are disturbed spirals (Bournaud et al. 
2007). This scenario seems to describe the observed mor¬ 
phologies of local LIRGs well. 

High-z ULIRGs differ from local objects with similar lu¬ 
minosities for several reasons. First, the incidence of merg¬ 
ers in high -z ULIRGs is lower than locally, with only 30- 
40% of the z > 1 ULIRGs showing merger morphologies 
(e.g., Elbaz et al.|2007 Kartaltepe et al.|2010 ). In addition, 


1 









































































Pereira-Santaella et al.: Star-formation histories of local LIRGs 


the mid-IR spectra of z ~ 2 ULIRGs differ from those of 
local U LIRGs because they are more similar to those local 
LIRGs (Farrah et ah, 2008 Rigby et ah, 2008). Therefore, 
the triggering mechanisms and the physical conditions of 
the SF in distant ULIRGs resemble those of local LIRGs. 
In this context, the detailed study of local LIRGs is needed 
to better understand their high -z counterparts. 

In this work we model the integrated spectral energy 
distribution (SED) of a sample of local IR bright galax¬ 
ies and derive fundamental physical parameters such as 
the stellar mass, SFR, star-formation history (SFH), av¬ 
erage extinction, etc. The analysis of the integrated emis¬ 
sion of local LIRGs is important for providing meaning¬ 
ful comparisons with the integrated SED used in high-z 
studies. Previous works of the integrated SED of nearb y 
galax i es focus on lower-luminosity galaxies (e.g ., Noll et ah] 


2009 Skibba et al. 2011 


Brown et al. 


2014) or high* 


luminosity U/LIRGs (e.g., jda Cnnha et al. 2010 U et al. 


2012 ), which shows that studying intermediate luminosity 
objects is needed. 

The paper is organized as follows. In Section |2 the sam¬ 
ple of local LIRGs is presented. Sections [3] and [4 describe 
the data reduction and the models used for the SED fitting, 
respectively. In Section [5] we discuss the SFH of this sam¬ 
ple and the age effects on the Ha to IR luminosity ratio. 
Finally, the main conclusions are presented in Section [6] 

Throughout this paper we assume the following cosmol¬ 
ogy: Hq = 70kms -1 Mpc -1 , U m = 0.3, and = 0.7. 


2. Sample 

We drew a volume-limited sample of local LIRGs from the 
IRAS Revised Bright Galaxy Sample (RBGS;|Sanders et al. 
2003). Our selectio n cr iteria are similar to those used by 


Alonso-Herrero et al. (2006|): Vh e i = 2750 — 5200 km s 
and Galactic latitude \b\ > 5, but we slightly decreased 
the minimum Lir down to log Lir/I/ 0 = 11.0. There are 
59 sources in the RBGS that fulfill these criteria, and 37 
of them are observable from the Roque de los Muchachos 
Obsevatory (Dec. > —16°). 

We obtained good quality g, r, and narrow-band Ha 
images for 25 of th em u sing the Nordic Optical Telescope 
(NOT; see Section 3.1), and anot h er fo ur had inte grated 
Ha flux measurements ( Moustakas fc Kennicutt[ |2006| ) and 
Sloan Digital Sky Survey (SPSS; Xihara et al.||201ip g and 
r images available. Therefore, our sample includes 29 (78%) 
of the parent sample northern LIRGs. Most of the missing 
galaxies belong to the lower end of the Lir distribution 
(log Lir/L 0 11.0), which is already well represented in 
our sample (see Table [l]). 

In addition, we observed six nearby companions of the 
RBGS objects (namely, NGC 876, UGC 03405, NGC 2389, 
MCG +02-20-002, NGC 6921, and NGC 7769), with 
logZ/iR/Z/Q= 10.2 — 10.7. They are located between 1 — 6' 
(20 — 100 kpc) away from the main RBGS galaxy and there¬ 
fore might contribute to the measured IRAS fluxes. We 
were also able to resolve three of the RBGS targets into two 
subcomponents (CGCG 468-002 NED01/02, NGC 7752/3, 
and NGC 7770/1). Our sample contains the 38 sources 
listed in Table 111 

The IR luminosity range is 10 10,2 -10 11,8 L©, with a 
mean and median luminosity of 10 110 L©. Except for one 
galaxy, CGCG 468-002 NED01, the IR luminosity is domi¬ 
nated by SF, and the bolometric AGN contribution is small, 


less than 5% for most galaxies and u p to 1 2% in a few of 
them (see Table [l] and Alonso-Herrero et al. 2012a). 

According to their nuclear activity classification, our 
sample includes 11 Hu galaxies, 13 composite, one LINER, 
five Seyfert galaxies, three objects without a clear opti¬ 
cal classificatiorQ and one galaxy with no available optical 
spectrum. For nine of them we determined their nuclear 
classification and [N n] 6584 A/Ha ratio using archival opti¬ 
cal spectroscopy (see Appendix |A|. The fraction of galaxies 
of each type is similar to what is expected for galaxies with 
IR luminosities in the range covered by our sample (Yuan 


et aU 2010). 


3. Observations and data reduction 

In this section we describe the reduc tion and analysis of 
our new optical obser vatio ns (Section |3.1| ) along with the 
archival data (Section |3.2| ). The reduced images for each 
galaxy are shown in Figure [l] and Appendix 0 


3.1. Optical imaging 

We obtained broad- and narrow-band imaging of 34 
IR bright galaxies using the Andalucia Faint Object 
Spectrograph and Camera (ALFOSC) on the 2.6 m NOT 
at the Roque de los Muchachos Obsevatory during three 
observing runs between May and December 2011 (see Table 
[2]) as part of the programs 115-NOT11/11A and 28- 
NOT2/11B. We used the broad-band SDSS g and r fil¬ 
ters (#120 and #110), and for the narrow-band images 
we used the filters #50 and #68 from the NOT filters 
set, and #65 and #66 from the Isaac Newton Group of 
Telescopes (ING). These narrow filters have a full-width 
half-maximum (FWHM) of about 5 nm and A c between 
663 and 665nm (see Tablej3|. For three objects (Arp 299, 
NGC 5936, and NGC 5990)not observed by us with the g 
filter we used their available SDSS g image (Aihara et al., 
2011). In addition, the r and g images of SDSS were used 
for IC 860, NGC 5653, Zw 049-057, and NGC 6052. 

We selected the narrow-band filter for each galaxy ac¬ 
cording to its redshift. The plate scale of ALFOSC is 
0.19"pixel -1 , and its field of view (FoV) is 7'x7', large 
enough to cover the emission of these LIRGs with a sin¬ 
gle pointing. The atmospheric conditions were photometric, 
and the seeing varied between O'/6 and l'/5 with a median 
seeing of 079. 

The integration times were 1300 and 800 s for the g and 
r filters, respectively, and 3000 s for the narrow-band filters. 
Each integration was divided into three to five dithered 
exposures that were later combined to correct for cosmic 
ray hits and bad pixels of the detector. 

For the data reduction we first subtracted the bias level 
using the overscan region to scale the master-bias. During 
the April run, the bias showed noticeable variability be¬ 
tween exposures, so we used the overscan region to estimate 
the bias level of each row in each exposure. Then the result¬ 
ing images were divided by the sky flat of the corresponding 
filter. In addition, bad pixels identified in the flat field im¬ 
ages were masked. The sky extinction was determined by 
observing spectrophotometric standard stars from the ING 
catalog (ING Technical Note 100) and standard SDSS stars 


1 Two or more lines of the 
not detected in their optica 


Baldwin et al. (1981) diagrams were 


spectra. 
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Table 1. The sample 


Name 

RA 

(J2000.0) 

Dec. 

(J2000.0) 

cz a 

(kms -1 ) 

D L b 

(Mpc) 

Spectral 

elass c 

[Nll]/Ha d 

Ref. 6 

log L\\i } 
(£©) 

Tagn/Tir^ 

NGC 23 

00 09 53.4 

+25 55 27 

4478 

64.7 

composite 

0.57 

1 

11.0 

0.02 

MCG +12-02-001 

00 54 04.0 

+73 05 05 

4722 

68.3 

H ii 

0.42 

1 

11.4 

<0.05 

NGC 876 

02 17 53.2 

+14 31 19 

3916 

56.5 

composite 

0.57 

2 

10.4 


NGC 877 

02 17 59.7 

+14 32 39 

3963 

57.2 

composite 

0.52 

3 

11.1 


UGC 01845 

02 24 08.0 

+47 58 11 

4600 

66.5 

composite 

0.72 

1 

11.1 

<0.05 

NGC 992 

02 37 25.5 

+21 06 02 

4065 

58.7 

H II 

0.43 

4 

11.0 


UGC 02982 

04 12 22.6 

+05 32 50 

5354 

77.6 

H II 

0.43 

5 

11.2 

<0.05 

NGC 1614 

04 34 00.0 

-08 34 45 

4778* 

69.1 

composite 

0.60 

5 

11.7 

<0.05 

CGCG 468-002 NED01 

05 08 19.7 

+17 21 48 

5267 

76.3 

Syl.9 

0.98 

2 

10.6 

0.40 

CGCG 468-002 NED02 

05 08 21.2 

+17 22 08 

4951 

71.6 

composite 

0.51 

2 

11.0 

<0.05 

NGC 1961 

05 42 04.7 

+69 22 43 

3908 

56.4 

LINER 

1.96 

6 

11.1 


UGC 03351 

05 45 47.9 

+58 42 03 

4433 

64.1 

Sy2 

1.34 

7 

11.2 

0.02 

UGC 03405 

06 13 57.6 

+80 28 35 

3799 

54.8 

composite? 

0.75 

2 

10.3 

<0.05 

UGC 03410 

06 14 29.6 

+80 27 00 

3871 

55.9 

H II 

0.45 

2 

10.9 

<0.05 

NGC 2388 

07 28 53.5 

+33 49 09 

4078 

58.9 

H II 

0.56 

1 

11.2 

<0.05 

NGC 2389 

07 29 04.7 

+33 51 40 

3956* 

57.1 

composite 

0.53 

2 

10.5 


MCG +02-20-002 

07 35 41.5 

+11 36 44 

5100* 

73.8 




10.2 


MCG +02-20-003 

07 35 43.4 

+11 42 35 

4907 

71.0 

composite 

0.45 

1 

11.0 

<0.05 

NGC 3110 

10 04 02.1 

-06 28 30 

5013 

72.6 

H II 

0.42 

5 

11.3 

<0.05 

NGC 3221 

10 22 20.0 

+21 34 10 

3959 

57.1 

H II 

0.40 

8 

11.0 


Arp 299 

11 28 31.0 

+58 33 41 

3056 

44.0 

Sy2 

0.35 

9 

11.8 

0.04 

MCG -02-33-098 

13 02 19.6 

-15 46 04 

4713 

68.2 

composite 

0.33 

5 

10.9 

0.16 

IC 860 

13 15 03.6 

+24 37 08 

3858 

55.7 

no 

7.8 

1 

11.0 

<0.05 

NGC 5653 

14 30 10.5 

+31 12 55 

3512 

50.6 

H II 

0.38 

3 

11.0 

0.01 

Zw 049-057 

15 13 13.1 

+07 13 32 

3858 

55.7 

H II 

0.46 

8 

11.0 

<0.05 

NGC 5936 

15 30 00.9 

+12 59 21 

3989 

57.6 

H II 

0.48 

1 

11.0 

0.03 

NGC 5990 

15 46 16.4 

+02 24 56 

3793 

54.7 

Sy2 

0.74 

10 

11.0 

0.05 

NGC 6052 

16 05 13.0 

+20 32 33 

4739 

68.5 

Hu 

0.23 

3 

10.9 


NGC 6701 

18 43 12.5 

+60 39 12 

3895 

56.2 

composite 

0.67 

1 

11.0 

<0.05 

NGC 6921 

20 28 28.9 

+25 43 24 

4329 

62.5 

AGN 


11,12 

10.2 

<0.05 

MCG +04-48-002 

20 28 35.1 

+25 44 00 

4198 

60.6 

Hn/AGN 

0.42 

11 

11.0 

0.06 

NGC 7591 

23 18 16.3 

+06 35 09 

4907 

71.0 

composite 

0.85 

1 

11.0 

<0.05 

NGC 7679 

23 28 46.7 

+03 30 41 

5161 

74.7 

Sy2/Syl 

0.59 

5 

11.1 

0.18 

NGC 7752 

23 46 58.5 

+29 27 32 

4943 

71.5 

composite 

0.33 

2 

10.7 


NGC 7753 

23 47 04.8 

+29 29 01 

5201 

75.3 

composite? 

1.02 

2 

10.9 


NGC 7769 

23 51 04.0 

+20 09 02 

4157 

60.0 

composite 

0.56 

2 

10.9 

0.10 

NGC 7770 

23 51 22.6 

+20 05 49 

4127 

59.6 

H II 

0.40 

13 

10.4 

0.17 

NGC 7771 

23 51 24.9 

+20 06 43 

4276 

61.8 

H II 

0.55 

1 

11.3 

0.02 


Notes. Heliocentric velocity from the Spitzer /IRS high-resolution spectra (Alonso-Herrero et al. 2012a). ^ Luminosity distance 
es timated from the red shift. ^ Classification of the nuclear activity from optical spectroscopy based on the classification scheme 
of Kewley et al. (2006). ^ d,e ^ Nuclear [Nil]6584 A/Ha ratio and reference. ^ Logarithm of the total 4-1000/im IR luminosity in 


solar units calculated in this pap er. The AGN toms IR lumino sity is included. ^ Ratio between the bolometric AGN luminosity 
and the total IR luminosity from Alonso-Herrero et al. (2012a). ^ Heliocentric velocities from NED. 


References. (1) Alonso-Herrero et al. 


(5) |Veill eux et al.| (|1995| ); (6) |Ho et al 

Kewley et al.| ( |2001| ); (11) Masetti et a 


A l); (3) |Moustakas &; Kennicutt ( 2006|); (4) Keel| ( 1984| ) ; 

I Aihara e t al.||20 1 1[ (9) |Garcfa-]V arm et al. ( |2006 ); (10) 


(2009); (2) This work ( Appen dix , . 

(1997); (7) |Baan et al. (|1998|); (8) Aihara et al.||2 01 1| 

. ( |2006| ); (12) |Tueller et al.| ( |2008 ); (13) Alonso-Herrero et al.| ( |2012b| ). 


Smith et al. 2002) at different air masses between 1 and 
Tdesides this, we used these standard stars to calculate 
the photometric AB zero points for the g and r filters. 


We combined individual exposures after subtracting the 
sky background emission and aligning them to a common 
reference image using stars in the FoV. The absolute as¬ 
trometry of each imag e was deter mined using the Guide 
Star Catalog (GSC2.3; Lasker et al. 20081). About 10-40 
stars were used for each image, and the estimated la astro- 
metric uncertainty is O'/1-0'/2. We checked our absolute flux 
calibration for those fields that were also observed by the 
SDSS (15 out of 28). We found good agreement between 
our measurements and those reported in the Eighth Data 


Release of the SDSS (Aihara et al. 2011) with differences 
around 0.01-0.05 mag for objects brighter than 19 mag. 

We used the Fitzpatrick (1999) extinction law and 
the Galactic color excess E(B — V ) (Table [5| from the 
NASA/IPAC extragalactic database (NED) to correct the 
observed optical fluxes for Galactic extinction. 


3.1.1. Ha narrow-band imaging calibration 

To obtain the Ha+[Nn] emission images, we first sub¬ 
tracted the continuum using the r band image. It was scaled 
using the relative fluxes of 10-30 field stars observed in both 
images using 4" apertures. These scaling factors are in good 
agreement with the theoretical factor that can be calculated 
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NGC 23 



NOT Ha+[NII] 


2MASSK 


50 

25 

0 

-25 

-50 


GALEX FUV JL GALEX NUV 


MIPS 24um 


PACS 70).im 


NOT g 


IRAC 3.6gm 


PACS 160|xm 


NOTr 


IRAC 8.0|i,m 


SPIRE 250|xm 


50 25 0 -25 -50 

RA (arcsec) 


Fig. 1. Subset of the images used to construct the integrated SED of NGC 23. The new NOT optical observations 
of the g , r, and Ha+[Nn] bands are shown. The rest of the images were obtained from the public archives of their 
respective observatories (see Section [3] for details). Some images used in the SED are not displayed here because they 
have morphologies similar to those of the presented images (2MASS J and H , IRAC 4.5 and 5.8 /im, and PACS 100 /im) 
or because they have a very low angular resolution (SPIRE 350 and 500/im). All images are shown in a logarithm scale. 
North is up and east to the left. The white bar in the H<a+[N n] panel represents lOkpc at the distance of the object. 


Table 2. Log of the NOT/ALFOSC observations 


Date 

Seeing 

(") 

Targets 

11 May 2011 

0.9-1.4 

Arp 299, MCG -02-33-098, NGC 5990 

12 May 2011 

0.9 

NGC 5936 

18 Sept. 2011 

0.8 

(NGC 6921, MCG +04-48-002), NGC 7591 

19 Sept. 2011 

0.7-1.1 

MCG +12-02-001, NGC 6701, NGC 7769, 
(NGC 7770, NGC 7771) 

29 Nov. 2011 

0.7-1.4 

NGC 23, (NGC 876, NGC 877), UGC 02982, 
NGC 1614, NGC 2388, NGC 2389 

30 Nov. 2011 

0.6-1.5 

UGC 01845, UGC 03351, MCG +02-20-002, 
MCG +02-20-003, NGC 3110, NGC 7679 

1 Dec. 2011 

1.2-1.5 

NGC 992, (CGCG 468-002 NED01/02), 
NGC 1961, (UGC 03405, UGC 03410), 

NGC 3221, (NGC 7752, NGC 7753) 


Notes. Galaxies observed with the same telescope pointing are grouped with parentheses. 


using the transmission curves of the filters. Some of the 
nights we obtained observations of standard stars with all 
the narrow filters, which confirms that the derived scaling 
factors are accurate within ~15%. 

To calculate the conversion factor from counts s -1 in the 
narrow filter to the corresponding Ha flux in erg cm -2 s -1 


units, we first created a synthetic spectrum including only 
the Ha line and the 6548 and 6584 A [N n] transitions at the 
redshift of each source. For the nuclear regions (^ 3kpc) 
where the [Nil ]/Ha ratio would be more uncertain be¬ 
cause of the nuclear activity, we used the ratios derived 
from the nuclear spectrum of each target (Table [l]). For the 
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Table 3. Characteristics of the NOT/ALFOSC filters 


Filter # 

Ac 

(nm) 

FWHM 

(nm) 

50 

665.3 

5.5 

68 

664.0 

4.0 

110 r 

618.0 

148 

120 g 

480.0 

145 

65 (ING) 

662.6 

4.4 

66 (ING) 

664.5 

5.0 


extra-nuclear emission we assumed that [Nnl6584A/ Ha = 


0.3, typical of Hu regions (Kennicutt & Kent, 1983). The 


synthetic spectrum was convolved with the transmission 
curve of the narrow filter and the result was converted into 
physical units using the known input Ha flux and the re¬ 
lation between the narrow- and broad-band r calibration 
(see previous section). The lOcr sensitivity of the images is 
^lO -16 erg cm -2 s _1 , which corresponds to an Ha luminos¬ 
ity of 10 37 ergs -1 at the median distance of our galaxies. 


archive to make the photomet ric me asurements using the 
apertures described in Section 3.1.2 (Table [4| taking into 
account that the angular resolution of the GALEX images 
is 4-6". As for the optical images, we masked stars inside 
of the apertures since in some cases they were bright, par¬ 
ticularly in the near-UV band. To estimate the background 
within the apertures, we used the sky background images 
provided by the GALEX pipeline. Finally to convert from 
count rates to physical units, we used the conversion factors 
given in the GALEX Observer’s Guide. 

For one of the galaxies without any GALEX data 
(NGC 1614), we took the near-UV flux from imaging ob¬ 
tained by the optical monitor (OM) onboard XMM- Newton 
usin g the UVW2 (2120A) filter (Pereira-Santaella et al. 


2011 ). 


We corrected the UV fluxes (Table [5| for Galactic ex¬ 
tinction using th e sam e method as we used for the optical 
data (see Section 3.1). 


3.2.2. 2MASS data 


3.1.2. Integrated photometry 

We defined the apertures to obtain the integrated emission 
using the NOT r images. In these images, we considered all 
the pixels above a surface brightness of 23 mag arcsec -1 and 
fitted an ellipse to them. The resulting elliptical apertures 
for each galaxy are listed in Table [4j They encompass most 
of the r emission for all the sources, although in a few cases 
some faint H II regions at large galactocentric distances (r > 
15kpc) lie outside of them. The diameters of the apertures 
range between 30 and 200" with a median diameter of 80", 
and are equivalent to ^20-25 kpc on average. 

To perform the photometry in the optical images (r, g, 
and Ha) we integrated all the emission in the calculated 
apertures (Table ffl) after masking Galactic stars that lie 
within them. The integrated fluxes are shown in Table |5j 

In this paper we use images with a wide range of angular 
resolutions, from ^1 to 35" (see next section). Therefore, to 
adapt the calculated elliptical apertures we convolved them 
using a Gaussian with a FWHM equal to the difference 
between the desired angular resolution and the r resolution 
(FWHM^l") subtracted in quadrature. 

3.2. Ancillary archival data 


To construct the SEDs of our sources, we looked for obser¬ 
vations carried out at different wavelengths from the UV 
to th e far-IR. In particu lar, we used image s from GALEX 


(UV; Martin et al.|2005 

), 2MASS (near-IR; 

Skrutskie et al. 

2006), Spitzer (mid-lR; Werner et al. 

2004 

), and Herschel 

(far-IR; Pilbratt et al. 2010p, all publicly available in their 


archives. The reduction and photometry of these archival 
data is described below. 


3.2.1. GALEX UV data 

In the GALEX archive we found far-UV (1516A) and/or 
near-UV (2367A) observations for 36 out of 38 galaxies in 
our sample. Most of them belong to the all-sky imaging sur¬ 
vey (AIS). The rest are part of the medium imaging survey 
(MIS), the nearby galaxies survey (NGS), or guest investi¬ 
gator programs. We used the images downloaded from the 


We retrieved the J (1.2/im), H (1.7/im), and K s (2.2/xm) 
near-IR images of our galaxies from the Two Micron All Sky 
Survey (2MASS). Mos t of them were from the 2MASS ex¬ 


tended source catalog (Jarrett et al., 2000), although a few 


objects were part of the 2MASS large galaxy atlas (Jarrett 


et al. 2003). These downloaded images were already - flux 


calibrated with an angular resolution of ~2". 

The photometry on the images was done considering the 
same elliptical apertures (Table®) used for the other bands. 
In general, there is good agreement between our integrated 
measurements (Table [5]) and those reported in the 2MASS 
catalogs. 


3.2.3. Spitzer IRAC and MIPS imaging 

In the Spitzer archive we found imaging of our galax ies fo r 
the four IR AC bands (at 3.6, 4.5, 5.8, and 8.0 pm; Fazio 
et al. [20044 and the 24 pm MIPS band ( Rieke et al.||2004p . 

We retrieved the basic calibrated data (BCD) for these 
observations from the Spitzer archive. The BCD process¬ 
ing includes several corrections (e.g., flat field, linearization, 
and dark subtraction) and flux calibration based on stan¬ 
dard stars. We combined these BCD images into mosaics 
using the version 18 of the MOsaicker and Point source 
Extractor (MOPEX) software provided by the SSC using 
the standard parameters (see the MOPEX User’s Guide 
for details on the data reduction). The FWHM of the point 
spread function (PSF) of these images vary between T/7 
and 2'/0 for the IRAC bands and it is 5'/9 for the 24 pm 
MIPS images. 

To measure the integrated emission, we used the aper¬ 
tures listed in Table [4] and corrected for the lower angu¬ 
lar resolution of the Spitzer images. For the IRAC images 
we applied the extended source aperture correction or the 
point source aperture correction for very compact objects 
(see the IRAC Instrument Handbook). These corrections 
are about 20% of the measured flux. In the MIPS 24 pm 
images, some galaxies are very compact, so we applied the 
point source correction described in the MIPS Instrument 
Handbook. 

MCG +02-20-002/3 were not observed with MIPS, 
therefore we consider the IRAS 25 pm flux from Surace 
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et al. (2004) as an upper limit because of the limited IRAS 


angular resolution (1' at 25 fi m) other sources might con¬ 
tribute to the IRAS flux. The measured integrated fluxes 
are listed in Tables [5] and |6l 


3.2.4. Herschel PACS and SPIRE imaging 


Far-IR imaging of our galaxies taken with He rschel /PACS 
(70, 100, and 160/xm; Poglitsch et al. |2010| ) and SPIRE 
(250, 350, and 500 /xm; Griffin et al.||2010p were available 
in the Herschel archive. Most of them were part of the 
program “Herschel-GOALS: PACS and SPIRE Imaging of 
a Complete Sample of Local LIRGs” (OTl_dsanders_l, PI: 
D. Sanders). 

To produce the images from the downloaded raw data, 
we first used the Herschel interactive pipeline environment 
software (HIPE) version 11 to create the flux calibrated 
timelines for each bolometer of the detectors. HIPE also 
attaches the pointing informatio n to the timelines. Then, 
using Scanamorphos version 22 (Roussel, 2012), we com¬ 
bined and projected these timelines in a spatial grid and 
obtained the final images. For the three PACS bands, the 
FWHMs of the PSF are ~6", 7", and 11", and for the 
SPIRE bands they are 18", 24", and 35", respectively. 

We performed the photometry on the PACS and SPIRE 
images using the apertures of Table [4] convolved with a 
Gaussian to account for the lower angular resolution of 
them. In some cases the galaxies were point-like at the 
Herschel resolution and we applied the point-source aper¬ 
ture corrections recommended in the PACS and SPIRE ob¬ 
server manuals. It was not possible to resolve the two galax¬ 
ies of CGCG 468-002 in the SPIRE 500 /xm image, so we 
used the measured flux as an upper limit to the emission 
of each component of the system. For NGC 7769 no PACS 
imag es were available, so we used the IRAS 100 /xm flux 
from Surace et al. (2004) as an upper limit estimate. In 
Table [6] we list the far-IR fluxes. 


4. SED modeling 

To investigate the properties of these galaxies from their 
SED, we used a metho d based on the models and fitting 
procedures presented by da Cunha et al. (2008). They com¬ 
pute the model SEDs from the UV to the far-IR wave¬ 
lengths taking the energy balance between the absorbed 
UV-optical radiation and that emitted in the IR by dust 
into account. First they calculate the emission from stars 
that is attenuated according to an extinction law, and this 
absorbed energy is re-emitted in the IR distributed into 
several components (PAH bands, hot grains, and warm and 
cold dust). These models are used together with the MAG- 
phys code (|da Cunha et al. 2008) to calculate the likeli¬ 
hood distributions of the physical parameters adopting a 
Bayesian approach. 

This method produces goo d results when applied t o 


nearby star-forming galaxies (da Cunha et al. 2008) 


However, we tried to use it directly with the SED of our 
LIRGs and th e results were not alw ays satisfactory as al¬ 
ready noted by da Cunha et al. (2010). This is because mod¬ 
els with physical conditions typical of LIRGs (higher ex¬ 
tinction and dust temperatures than in normal star-forming 
galaxies) are not numerous in their set, so the obtained like¬ 
lihood distributions are not reliable. Therefore, to analyze 


Table 4. Photometry apertures 


Name 

a° 

(") 

b a 

(") 

P.A 6 

O 

2 a c 
(kpc) 

NGC 23 

40.1 

32.5 

150 

25.2 

MCG +12-02-001 

16.7 

8.2 

117 

11.0 

NGC 876 

49.4 

17.0 

28 

27.1 

NGC 877 

55.5 

42.5 

139 

30.8 

UGC 01845 

24.6 

12.2 

139 

15.9 

NGC 992 

25.4 

15.5 

13 

14.4 

UGC 02982 

33.9 

14.7 

108 

25.5 

NGC 1614 

29.5 

24.1 

163 

19.8 

CGCG 468-002 NED01 

18.3 

13.0 

94 

13.6 

CGCG 468-002 NED02 

27.2 

4.6 

147 

18.9 

NGC 1961 

106.0 

68.8 

91 

58.0 

UGC 03351 

41.4 

9.2 

166 

25.7 

UGC 03405 

40.6 

9.4 

128 

21.6 

UGC 03410 

48.6 

13.4 

120 

26.3 

NGC 2388 

26.5 

17.6 

69 

15.1 

NGC 2389 

64.4 

40.1 

74 

35.6 

MCG +02-20-002 

24.9 

18.0 

94 

17.8 

MCG +02-20-003 

22.9 

12.8 

144 

15.8 

NGC 3110 

38.8 

18.9 

171 

27.3 

NGC 3221 

96.8 

25.6 

167 

53.6 

Arp 299 

43.5 

41.2 

7 

18.6 

MCG -02-33-098 

46.4 

11.5 

62 

30.7 

IC 860 

29.1 

16.5 

18 

15.7 

NGC 5653 

46.1 

33.5 

112 

22.6 

Zw 049-057 

17.8 

11.1 

24 

9.6 

NGC 5936 

42.0 

39.3 

156 

23.4 

NGC 5990 

37.0 

23.9 

116 

19.6 

NGC 6052 

27.2 

20.7 

4 

18.1 

NGC 6701 

36.3 

33.6 

27 

19.8 

NGC 6921 

34.1 

12.0 

140 

20.7 

MCG +04-48-002 

24.9 

10.1 

67 

14.6 

NGC 7591 

48.6 

24.1 

143 

33.4 

NGC 7679 

26.5 

19.4 

79 

19.2 

NGC 7752 

20.6 

9.7 

98 

14.3 

NGC 7753 

56.6 

43.7 

76 

41.3 

NGC 7769 

43.2 

35.5 

123 

25.2 

NGC 7770 

20.5 

15.6 

12 

11.9 

NGC 7771 

53.1 

21.0 

71 

31.8 


Notes. Elliptical apertures used to measure the integrated emis¬ 
sions based on the 23magarcsec -1 isophotes. ^ Semi-major 
and -minor axes of the elliptical aperture. ^ Position angle 
measured counter-clockwise from the North axis. Physical 
size of the major axis of the aperture at the assumed distance 
(see Table [I]). 


our data better and to include the H<a emission in the fit, we 
decided to generate a new set of models for the stellar and 
IR emissions and to modify the original magphys code. In 
the following we describe our models and modifications. 


4.1. Stellar emission 


We modeled the stellar emission by combining the emission 
of single stellar population bursts of different ages weighted 
by the stellar mass of each age. The input ste l lar spectra are 


from the POPSTARp] library (Molla et al. 2009 Martin- 
Manjon et al. 12010). Specifically, we used only those models 


created with the 

and solar metallicity. 


Kroupa (2001) initial mass function (IMF) 


2 http://www.fractal-es.com/PopStar 
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Table 5. Integrated UV, optical, and near-IR photometry 


Name 

E(B - V) a 
(mag) 




Integrated fluxes (mjy) 




H a b 
6563A 

GALEX 

NOT 


2MASS 


Spitzer / IRAC 

FUV 

1516A 

NUV 

2267A 

9 

4770A 

r 

6230A 

J 

1.23/rm 

H 

1.66/zm 

K 

2.16/im 

3.6/xm 

4.5gm 

NGC 23 

0.034 

13 

1.2 

2.8 

37 

67 

170 

210 

180 

100 

68 

MCG +12-02-001 

0.550 

14 

<1.80 

<2.94 

7.4 

14 

47 

58 

56 

51 

45 

NGC 876 

0.099 

1.0 

0.17 

0.30 

5.0 

9.8 

34 

45 

45 

26 

19 

NGC 877 

0.099 

24 

4.3 

7.5 

61 

99 

220 

260 

220 

130 

93 

UGC 01845 

0.190 

2.3 


<0.17 

6.1 

14 

55 

80 

78 

50 

37 

NGC 992 

0.130 

14 

1.5 

2.7 

18 

30 

76 

96 

84 

63 

46 

UGC 02982 

0.360 

7.7 

<0.57 

1.1 ± 0.2 

12 

22 

68 

86 

85 

62 

46 

NGC 1614 

0.140 

18 


2.9 C 

26 

41 

97 

120 

120 

93 

76 

CGCG 468-002 NED01 

0.310 

1.9 

<0.35 

0.47 

11 

19 

50 

63 

59 

42 

41 

CGCG 468-002 NED02 

0.310 

1.7 

<0.37 

0.70 

4.0 

6.6 

17 

22 

19 

12 

11 

NGC 1961 

0.110 

26 

4.9 

10 

100 

170 

470 

600 

500 

260 

170 

UGC 03351 

0.250 

3.0 



7.2 

16 

76 

110 

120 

77 

58 

UGC 03405 

0.085 

1.5 

0.15 ± 0.02 

0.31 

5.7 

12 

40 

53 

50 

30 

20 

UGC 03410 

0.085 

5.7 

0.40 

0.75 

11 

23 

86 

120 

110 

78 

57 

NGC 2388 

0.051 

4.0 

0.049 

0.21 

9.7 

22 

78 

100 

95 

64 

48 

NGC 2389 

0.051 

11 

3.5 

5.4 

23 

33 

60 

75 

59 + 6 

37 

26 

MCG +02-20-002 

0.028 

2.7 

0.95 

1.6 

8.5 

12 

22 

24 

21 

12 

8.3 

MCG +02-20-003 

0.028 

3.0 

0.42 

0.79 

6.5 

11 

25 

32 

30 

26 

33 

NGC 3110 

0.031 

14 

1.6 

3.3 

25 

40 

92 

110 

100 

75 

54 

NGC 3221 

0.021 

7.1 

0.90 

1.7 

24 

45 

140 

210 

180 

120 

83 

Arp 299 

0.015 

80 

9.2 

13 

59 d 

95 

210 

270 

270 

290 

370 

MCG -02-33-098 

0.053 

4.9 

0.60 

1.2 


19 

45 

55 

49 

33 

26 

IC 860 

0.012 

0.12 e 

<0.02 

0.19 

9.0 d 

16 d 

37 

48 

38 

19 

13 

NGC 5653 

0.013 

14 e 

1.6 

3.4 

33 d 

58 d 

130 

160 

140 

92 

66 

Zw 049-057 

0.035 

0.51 e 

<0.01 

0.057 

3.4 d 

7.2 d 

21 

29 

24 

14 

11 

NGC 5936 

0.034 

12 

2.3 

4.2 

31 d 

48 

110 

130 

110 

72 

51 

NGC 5990 

0.099 

14 

1.5 

3.4 

31 d 

54 

150 

180 

170 

130 

120 

NGC 6052 

0.067 

26 e 

4.5 

7.0 

22 d 

28 d 

39 

44 

36 

28 

21 

NGC 6701 

0.037 

11 


3.8 

32 

56 

130 

150 

130 

81 

57 

NGC 6921 

0.390 

1.4 


<0.89 

40 

76 

190 

240 

210 

96 

61 

MCG +04-48-002 

0.390 

6.1 


<1.04 

10 

20 

64 

85 

81 

60 

52 

NGC 7591 

0.091 

5.6 

0.53 

1.3 

17 

33 

89 

110 

97 

55 

39 

NGC 7679 

0.058 

16 

2.2 

4.3 

25 

36 

65 

73 

64 

52 

39 

NGC 7752 

0.088 

6.8 

0.93 

1.7 

7.6 

11 

23 

26 

24 

19 

13 

NGC 7753 

0.087 

10 

2.5 

4.6 

40 

68 

170 

220 

190 

84 

56 

NGC 7769 

0.066 

11 

3.9 

6.1 

42 

73 

170 

210 

170 

95 

63 

NGC 7770 

0.066 

6.9 

1.3 

2.0 

9.6 

16 

30 

37 

31 

25 

24 

NGC 7771 

0.066 

10 

0.87 

2.1 

36 

73 

230 

300 

270 

150 

110 


Notes. Statistical uncertainties are included only when they are larger than 10% of the flux. UV, optical, and Ho fluxes are 
corrected for Galactic extinction using the Fitzpatrick (1999) extinction law. Galactic color excess E(B — V ) from NED used 
for the Galactic ext inction correction. ^ The units of the observed Ho fluxes are 10 -13 erg cm -2 s -1 . ^ XMM-Newto n /OM UVW2 


(2120 A) fl u x fro m 


Kennicutt (2006). 


Pereira-Santaella et al. 


(2011). ^ Measured from SDSS images. ^ Integrated Ho fluxes from 


Moustakas & 


In this library there is a large number of models (106) 
for different stellar ages but many of them have very sim¬ 
ilar photometric colors. Because of this, we selected a set 
of four representative spectra of several age ranges (10- 
100 Myr, 100-500 Myr, 0.5-1.5Gyr, and >1.5 Gyr; stars 
younger than 10 Myr are considered below through a re¬ 
cent star-formation burst, see below). We chose these age 
ranges because the variation in the mass-to-luminosity ra¬ 
tio is lower than a factor of two, and they have similar 
photometric colors, so age variations within these ranges 
would be almost indistinguishable using photometric infor- 
matio n alone. We c ompared the POP STAR models with 
the Maraston (2005) models, which include a detailed treat- 
ment of the thermally pulsing asymptotic giant branch 
(TP-AGB). The main differences appear in the near-IR 
range for populations with 0.5-1.5 Gyr, although the dif¬ 
ferences in luminosity are small, a factor of 2-3. 


To combine the single stellar population models, we 
considered that a random fraction of the stellar mass was 
formed at a constant rate during the five age intervals men¬ 
tioned above. We also added to the SFH a recent burst of 
star-formation that began between 1 and 300 Myr ago (£sb) 
and continues until today with a constant intensity I$b be¬ 
tween 0.03 and 10 000 times the average previous SFR. 


Then, we extinguished the com bined stellar sp e ctrum 
following the prescription given by Chariot & Fall (2000) 
and used by |da Cunha et al. (2008) in the original MAG- 
phys code. That is, we assumed that stars younger than 
10 Myr are still embedded in their birth clouds (BC) and 
have higher extinctions than older stars, which are only af¬ 
fected by the interstellar medium (ISM) extinction. Thus, 
for a given wavelength, the total effective absorption, ry, 


is r; 


BC 


I MSM 
T r A 


for stars younger than 10 Myr and rj( 


ISM 


for older stars. For the wavelength dependence of ry, we 


used the power-law dependence assumed by Chariot fc Fall 
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Table 6. Integrated mid- and far-IR photometry 


Name 

Spitzer/ IRAC 

Spitzer/ MIPS 


Herschel /PACS 


Herschel/SPIRE 

5.8 /im 

8.0 /im 

24 /im 

70 /xm 100 /im 

160 /im 

250 /im 

350 /im 

500 /im 

NGC 23 

0.17 

0.43 

0.86 

10 

14 

13 

5.2 

2.1 

0.63 

MCG +12-02-001 

0.20 

0.58 

2.9 

25 

29 

22 

7.5 

2.8 

0.96 

NGC 876 

0.045 

0.10 

0.17 

3.3 

5.0 

5.5 

2.6 

1.1 

0.39 

NGC 877 

0.29 

0.77 

0.88 

12 

23 

28 

14 

5.8 

1.9 

UGC 01845 

0.12 

0.34 

0.84 

13 

17 

15 

4.8 

1.8 

0.56 

NGC 992 

0.19 

0.55 

1.0 

12 

16 

15 

5.4 

2.1 

0.66 

UGC 02982 

0.19 

0.56 

0.60 

9.9 

16 

16 

7.2 

3.1 

1.2 

NGC 1614 

0.30 

0.88 

5.7 

36 

36 

23 

6.5 

2.3 

0.65 

CGCG 468-002 NED01 

0.056 

0.088 

0.33 

2.4 

2.9 

2.5 

0.82 

0.30 

<0.19 a 

CGCG 468-002 NED02 

0.033 

0.086 

0.60 

9.2 

9.6 

6.8 

2.0 

0.73 

<0.19 a 

NGC 1961 

0.32 

0.73 

0.71 

11 

24 

36 

18 

8.1 

3.0 

UGC 03351 

0.19 

0.52 

0.66 

18 

31 

32 

14 

5.5 

1.6 

UGC 03405 

0.054 

0.14 

0.12 

2.1 

4.3 

5.6 

2.9 

1.3 

0.45 

UGC 03410 

0.19 

0.51 

0.62 

9.9 

17 

20 

8.8 

3.6 

1.2 

NGC 2388 

0.15 

0.42 

1.6 

20 

25 

21 

7.4 

3.0 

0.97 

NGC 2389 

0.079 

0.20 

0.27 

3.2 

6.3 

5.2 

2.3 

1.1 

0.50 

MCG +02-20-002 

0.021 

0.057 

C0.50 6 

0.85 

1.6 

2.1 

0.96 

0.43 

0.16 

MCG +02-20-003 

0.11 

0.22 

C0.73 6 

10 

13 

10 

3.5 

1.4 

0.45 

NGC 3110 

0.21 

0.60 

0.91 

13 

20 

20 

7.9 

3.1 

0.94 

NGC 3221 

0.24 

0.64 

0.72 

10 

19 

25 

12 

5.1 

1.7 

Arp 299 

1.1 

2.4 

23 c 

130 

120 

72 

22 

7.6 

2.3 

MCG -02-33-098 

0.076 

0.21 

1.1 

7.2 

8.4 

6.5 

2.6 

1.1 

0.34 

IC 860 

0.018 

0.038 

0.87 

20 

18 

11 

3.6 

1.5 

0.47 

NGC 5653 

0.22 

0.62 

1.1 

14 

21 

21 

8.1 

3.1 

1.0 

Zw 049-057 

0.028 

0.071 

0.54 

28 

32 

24 

7.9 

3.0 

0.98 

NGC 5936 

0.17 

0.48 

1.0 

11 

16 

16 

6.7 

2.6 

0.87 

NGC 5990 

0.23 

0.48 

1.2 

12 

16 

15 

6.3 

2.4 

0.84 

NGC 6052 

0.075 

0.20 

0.66 

7.4 

9.9 

8.9 

3.3 

1.3 

0.43 

NGC 6701 

0.17 

0.45 

0.98 

12 

19 

18 

7.4 

2.9 

0.97 

NGC 6921 

0.060 

0.077 

0.086 

2.2 

3.1 

2.3 

0.79 

0.30 

0.10 ± 0.03 

MCG +04-48-002 

0.17 

0.46 

0.68 

10 

14 

12 

4.1 

1.5 

0.49 

NGC 7591 

0.10 

0.28 

0.65 

8.9 

13 

13 

6.4 

2.7 

0.91 

NGC 7679 

0.14 

0.40 

0.81 

8.0 

10 

8.4 

3.1 

1.2 

0.42 

NGC 7752 

0.052 

0.14 

0.28 

3.4 

4.9 

4.5 

1.7 

0.70 

0.21 

NGC 7753 

0.11 

0.26 

0.34 

4.1 

8.6 

12 

6.5 

3.0 

1.1 

NGC 7769 

0.14 

0.35 

0.49 


<13.58 6 


4.9 

1.9 

0.61 

NGC 7770 

0.060 

0.14 

0.38 

2.9 

3.8 

3.2 

1.2 

0.44 ± 0.07 

0.14 ± 0.03 

NGC 7771 

0.26 

0.63 

1.2 

23 

38 

39 

17 

6.6 

2.1 


Notes. The units of the fluxes are Jy. Uncertainties are included only when they are larger than 10% of the flux. ^ The two 
components of CGCG 468-002 are not resolved in the 500 /im Herschel /SPIRE image, therefore we assume the inte grated flux of 
t he tw o components as an upper limit. ^ These upper limits are the IRAS 25 /im and 100 /im flu xes measured by 
(2004). The MIPS 24 /im image of Arp 299 is saturated so we took the I RAS 25 /im flu x from 
by the conversion factor between IRAS 25 /im and MIPS 24 /im fluxes given Calzetti et al. (2010). 


Sanders et al. 


Surace et al. 


(2003) corrected 


(2000). Similar to da Cunha et al. (20081), the input parame¬ 
ters regarding the extinction are r v and /i, where T+uis the 
total extinction affecting young stars 


and n = + M / t, 


We also computed the amount of energy that is absorbed 
by dust (Adust) and the fraction of this absorbed energy 
that is produced by stars older than lO Myr (J M ). These 
two parameters are used later in Section |4.4| to combine 
the stellar models with the IR emission models. 

The luminosities of the hydrogen recombination lines 
were calculated from the number of ionizing photons in the 
unextinguished stellar spectrum using the case B |Storey fc 


Hummer (1995) recombination coefficients. These emission 
lines are affected by the BC and/or the ISM extinctions as 
well. 


3 A v /r v = 2.5 log e ~ 1.086. Both A v and r v are used inter¬ 
changeably along this paper taking into account this factor. 


In total we generated 50 000 different models, which is 
enough to produce smooth likelihood distributions in the 
fits. 


4.2. IR emission 


For the IR emission we used a two-component model: dust 
thermal emission and PAH emission. For the former com¬ 
ponent we assumed that the temperature distribution of 
the dust mass follows a power-law, dM^ /dT oc T~ 7 with 
a low temperature cut-off, T m i n (see, e.g., Dale et aT||2001 


Kovacs et al. 2010|). We assumed that the dust emission for 


a given temperature is a graybody with fixed /3 = 2 and 
an abso rption coefficient k = 0.517 m 2 kg -1 at 240 /im (Li 


& Draine 2001). We estimated the fraction of the IR lumi¬ 
nosity produced by dust with T < 50 K (/^ R ) to separate 
the fraction of the IR luminosity that is produced in pho- 
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Fig. 2. Top panel : Best-fitting model for the SED of NGC 23 (constant SFH, see Section 5.1) in a solid black line, 
together with the observed data (red diamonds). The Ha luminosity and model prediction in L® units are plotted at 
A = 6563A (red diamond and solid black circle, respectively) multiplied by a factor of 10 3 . The solid color lines indicate 
the contribution of the different stellar populations with ages <10Myr (blue), 10-100Myr (dark green), 100-500Myr 
(light green), 0.5-1.5 Gyr (orange), and > 1.5 Gyr (pink). The blue and red dashed lines are the dust emission for dust 
tem peratures l owe r and higher than 50 K, respectively. The dotted green line is the AGN torus model derived by ||Alonso - 
|Herrer o et al. ( |2012a ). The residuals are shown in the lower part of the panel. Bottom panels: Likelihood distribution for 
several of the physical parameters (IR dust luminosity, SFR, stellar mass, ISM extinction /ir v , young stars extinction r v , 
ratio between the ISM and young stars extinctions /i, fraction of IR luminosity produced by cold dust with T < 50 K 
removing the AGN contribution /^, fraction of IR luminosity due to PAH emission, and dust mass). The next five panels 
show the logarithm of the percentage of the stellar mass for different stellar age ranges using the same color coding as in 
the top panel. The last panel represents the SFH. 


todissociation regions (PDR, U ~ 20Cp|) and that produced 4.3. AGN contribution 

20 ' 


in more diffuse regions (D raine fc Li[ 2007). 

The PAH emission consists of several emission bands in 
the mid-IR range between 3 and 20 /am. To obtain a “pure” 


PAH template we used the Smith et al. (2007) average 5 


35 /im mid-IR spectra of local star-forming galaxies after 
subtracting the underlying hot dust continuum using the 


pahfit code (Smith et al. 2007). We added the PAH fea¬ 
ture at 3.3 p m using a Drude profile wit h an intensity equa l 
to one third of the 6.2 pm PAH feature (Draine & Li 2007). 


random values for 7 , 
(<Zpah)- 


Most of the galaxies in our sam ple are part of the large r 
sample of LIRGs studied by Alonso-Herrero et al. ( 2012 a). 
In that work the AGN contribution to the mid-IR emis- 
sion was estimated decomposing their 5-38 pm Spitzer low- 
resolution sp ectra using star-for mation templates (Brandi 


et al. 


2006; 


The dust emission model and the PAH template were 
combined into the final model assuming that the PAH lu¬ 
minosity can represent between 1 and 40% of the total IR 
luminosity. We produced 20 000 IR emission models with 


els ( |Nen^o ~ 


iva 


Rieke et 
et al 


al. 


rpoosi). 

he IR is 


and PAH luminosity fractions 


2009) and clumpy torus mocC 
Although, in general, the AGN 
energy output in the IR is small compared to that of star- 
formation in these galaxies (see Table [lj, at certain wave¬ 
lengths the AGN contribution can be noticeable. For this 
reason we us ed the torus model fitted by |Alonso-Herrero 


et al. 
AGN 


2012a), when an AGN was detected, to subtract the 
R emission from our integrated measurements. 


4 Where U = 1 is the interstellar radiation field in the solar 
neighborhood. 


Only for CGCG 468-002 NED 01 and NGC 7770 does 
the AGN dominate the mid-IR emission. 
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Fig. 3. Same as Figure [ 2 ] but for Arp 299 (recent SF burst, see Section [ 5 TT] ) 


12 3 4 

log t (Myr) 


4.4. Bayesian parameter inference 


To determine the likelihood distributions of the physical 
parameters we modified the magphys code so we could 
make use of the models we constructed. First, stellar and IR 
models are combined to obtain the complete SED requiring 
that /^ R = fn dz 30%. In total, we find about 470 million 
combinations that fulfill this requirement. Then, the SED 
models are scaled to match the observed photometric fluxes 
and Ha emission. A probability ( e~ x / 2 ) is assigned to each 
model by considering upper limits when present in the SED 
(see Appendix H- Finally, the likelihood distributions of 
the parameters are derived from these probabilities (see da 
|Cunha et al.||2008| for details). 


4.5. Fitting results 


To fit the SED of our sample of LIRGs (Tables [ 5 ] and [6| 
we used the procedure described in the previous section. 
For most of these integrated measurements, the statisti¬ 
cal error is low (< 5%) so most of the uncertainty comes 
from systematic errors, such as the absolute flux calibra¬ 
tion, the uncertain aperture correction for semi-extended 
sources, and possible aperture mismatches. Therefore, we 
assumed a conservative 20% systematic error for all the 
photometric points in our SED added in quadrature to the, 
typically very low, statistical error. 


In Figures [2] and |3] we show the results of SED fitting 
for two of our galaxies: one with a constant SFH and the 
other with a strong burst of recent SF (see Section 5.1). For 


the rest of the sample, they are shown in Appendix [D] The 
fitted parameters are listed in Table [8] 

The best-fitting model for each galaxy is able to closely 
reproduce the observed SED for most of the galaxies except 
for IC 860 and Zw 049-057. Their optical g and r emis¬ 
sions are clearly underestimated by the best-fitting models 
(see Appendix |D|) . This can be caused by the presence of a 
deeply obscurea energy source (AGN or SF), which is only 
detected through its far-IR emission or an extinction be¬ 
havior that is more complex than the one assumed in the 
models. In either case, the derived parameters may be un¬ 
certain, so we d eci ded t o exclude them from the discussion 
of Sections 15.11 and 15.21 

5. Results and discussion 

5.1. Star-formation histories 

From the Bayesian fitting of the SEDs we derived likeli¬ 
hood distributions for the SFH of these galaxies based on 
the different SFH of the models. As explained in Section 
|4.1[ with the photometric information alone, we are able 
to distinguish stellar populations in limited age intervals. 
Because of this, we calculated the SFH in four intervals with 
the same duration in log scale (i.e., 0-10 Myr, 10-100 Myr, 
0.1-1 Gyr, 1-10 Gyr). These intervals are almost exactly the 
same as we used to construct the SED models. 

To classify these SFH we first tested whether they are 
compatible with a constant SFR (see Figure [4]). For 17 out 
of 36 galaxies we found that their SFH deviate less than 
3cr from a constant SFR during their histories (averaged 
over the time intervals specified above). The mean uncer- 
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Constant SFR SF burst Decaying SFR 



t (Gyr) t (Gyr) t (Gyr) 


Fig. 4. Star-formation histories derived from the SED fitting. In the left panel we show the SFH of those galaxies with 
a constant SFH in blue. In the central panel, we plot those SFH with a recent burst of SF (green). In the right panel we 
show those SFH with a decaying SFR (orange). The numbers represent the galaxy ID (see Table [9]). 



Fig. 5. SFR vs. stellar masses at three different age intervals calculated from the SFH. The color of the s ymbols is as 
i n Fig ure |4j The solid and dashed black lines are the M-S relation and uncertainty, respectively, derived by |Elbaz et al.| 
(2007) for SDSS galaxies at z^O. The horizontal dashed red line marks the approximate thr eshold SFR needed to reach 
a logLiR/Z /0 > 11 (excluding the AGN luminosity) based on our calibration (Section |5.2| ). Galaxies below the LIRG 
threshold are close companions of the sources identified as LIRGs using low angular resolution ( 2 ') IRAS data (see 
Section [ 2 | . 


tainty of the SFR in each interval is ^0.5 dex, so we are 
not able to detect variations lower than a factor of ^3 in 
the SFR. These 17 galaxies are mostly spirals, some of 
them with clearly disturbed morphologies (e.g., NGC 1961 
or NGC 5653). For two cases we found a decaying SFR, and 
the remaining 17 galaxies (50% of the sample) show a SFH 
with a SF burst in the 0-100 Myr intervals. We did not find 
significant differences between the SFR in the 0-10 Myr and 
10-100 Myr ranges for any of the galaxies, including those 
with a burst of SF (starbursts in the following) and those 
with nearly constant SFR. 


The SFH of the identified starbursts indicates that the 
current burst began on average between 30-300 Myr ago. 
The upper limit of this range is calculated by estimating 
how long it would be possible to sustain the burst intensity 
during the 0.1-1 Gyr interval and still detect a relative SF 
burst in the 10-100 Myr interval. Similarly, the lower limit 
is calculated by doing the same for the 1-10 Myr and 10- 
100 Myr int ervals. This burst dura tion is similar to the one 
derived by Marcillac et al. (2006), 40-260 Myr from the 
modeling of Balmer absorption lines and the 4000A break 


in z = 0.7 LIRGs. The intensity of the current burst is 
between 2 and 20 times with a median of 7.5 times the 
previous averaged SFR (Table [9]). 

Galaxies with the higher burst intensities (> 10) 
are those whose morphologies are highly disturbed in¬ 
dicating recent or ongoing interactions (MCG + 12 - 02 - 
001 , NGC 1614, CGCG 468-002 NED 02 , Arp 299, and 
NGC 6052). However, in our sample there is one object 
(MCG -02-33-098) that is also an interacting system with 
two nuclei separated by ^4kpc, but its burst intensity is 
only (2.5lo' g )- Combining the age of the burst with its 
intensity, we calculate the stellar mass produced by the 
current burst of SF with respect to the current stellar mass 
(Table [9|, which can be >10-30% in these mergers. For the 
rest of the galaxies the mass formed tends to be < 10 %, al¬ 
though there is a continuous distribution of formed mass 
fractions between the extreme bursts of mergers and the 
weaker SF bursts of other galaxies in our sample. 

In Figure[5]we plot the SFR averaged over the 0-10 Myr, 
10-100 Myr, and 0.1-1 Gyr intervals vs. the stellar mass to¬ 
day, 100 Myr ago, and 1 Gyr ago. From the left hand panel 
of this figure, it is clear that those galaxies classified as 
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starbursts lie above the SFR-M* main-sequence (M-S), as 
expected, while those with a relatively constant SFH are 
consistent with the M-S relation within 2 a. 

Figure [5] shows the luminosity threshold for a galaxy 
to be classified as LIRG (logLiR/L 0 >11). According to 
this luminosity criterion, there are 8 starbursts and 5 M-S 
galaxies with LIRG luminosities due to star-formation in 
our sample. The latter set of galaxies (NGC 23, NGC 877, 
NGC 1961, NGC 2388, and NGC 7771) are spirals classi¬ 
fied as LIRGs, although they do not have a particularly 
high specific SFR (sSFR = SFR/M* < 0.4 Gyr" 1 ). They 
lie within lcr in the M-S, and since they are relatively mas- 
sice (logM*/M© > 11.0), their expected SFR imply IR 
luminosities above or close to the LIRG luminosity thresh¬ 
old. 


The starburst LIRGs have enhanced SFR with respect 
to the M-S, with the most extreme cases (Arp 299, 
NGC 6052, MCG +12-02-001, and NGC 1614) being merg¬ 
ers with high specific sSFR > 1.0 Gyr 1 (see Table [8]). 
The rest of starbursts include spiral galaxies with differ¬ 
ent degrees of disturbed morphologies (e.g., NGC 6701 and 
NGC 5936). The AGN contribution to the IR lu minosities 
is sm all in our sample (Table [l] and Alonso-Herrero et al.| 
2012al, therefore the presence of an AGN does not affect 


their classification as LIRGs. 

Both the left hand and middle panels of Figure [5] 
show similar distributions of the galaxies in the M*-SFR 
plane. As indicated before, this is probably because the 
star-formation bursts have been active during the past 
^100 Myr. However, in the righthand panel (0.1-1 Gyr 
range; z ~ 0.075), the situation is completely different. 
Almost all the galaxies are in agreement with the M-S rela¬ 
tion, and only one or two objects, with log M*/M 0 > 11.1, 
would have been classified as LIRGs if they had been ob¬ 
served J 1 Gyr ago. 


5.1.1. Comparison with previous results for local U/LIRGs 


Rodriguez Zanrfn ~t al.| (|2009, |2010 ) used optical spec¬ 
troscopy to study the SFH of local ULIRGs and find that 
young stellar populations (age <100 Myr) usually dominate 
the stellar mass of these systems. Most of their ULIRGs 
are mergers in different evolutionary stages, so we compare 
their results with the five mergers in our sample with burst 
intensities > 10. In our merger LIRGs, the stellar mass 
formed during the ongoing burst is up to ~30%, slightly 
lower than in ULIRGs, although they are also compatible 
with burst stellar mass fractions as low as 5-10% (see Table 

0 - 

Rodriguez Zaurfn et al. (2009) argue that an old stellar 
population (> 2 Gyr) is not present in most of the optical 
spectra of ULIRGs, which is consistent with the best-fitting 
models for our merger LIRGs (Figure and Appendix Dj) 
where the optical light is dominated by young stars (10- 
100 Myr). 

Similar results regarding the stellar m asses of ULIRGs 
were obtained by da Cunha et al. (2010) using magphys. 
In that study they also found that the median dust mass 
in their ULIRGs is 1O 8,6 M 0 , which is almost a factor of 
10 higher than the dust masses derived for these LIRGs. 
However, the dust mass determination in ULIRGs is uncer¬ 
tain because it is not easy to constrain the cold dust temper¬ 
ature, which contributes sign ificantly to the dus t mass but 
not to the IR luminosity (da Cunha et al., 2010). Actually, 


Table 9. Star-formation burst properties 


Name 

Intensity a 

Age b 

Mass c 



(Myr) 

(%) 

MCG +12-02-001 

13+ s 

40-200 

3-16 

NGC 992 

6+ 5 

20-440 

1-15 

UGC 02982 

o 9+4.9 

z - z -o.i 

40-210 

1-3 

NGC 1614 

1( +0 3 7 

50-370 

4-27 

CGCG 468-002 NED 02 

12 > 

40-290 

4-25 

NGC 2389 

o i +4.0 

50-310 

2-8 

MCG + 02 - 20-002 

9 90+O.47 
Z - ZO -0.83 

30-190 

1-4 

MCG +02-20-003 

7.0+J 

40-170 

3-10 

NGC 3110 

3-Otii 

30-200 

1-6 

Arp 299 

17+ 4 

70-260 

9-29 

MCG -02-33-098 

9 k+5.4 
Z -°- 0.8 

40-140 

1-3 

NGC 5936 

2 . 1 ±?!o 

30-110 

1-3 

NGC 6052 

18.8+2 

60-420 

12-60 

NGC 6701 

1 - 8+8 

30-40 

1-1 

NGC 7679 

8-4+t 

60-290 

6-22 

NGC 7752 

7.9+;! 

30-210 

2-13 

NGC 7770 

4 +i 0 

20-260 

1-7 


Notes. Galaxies with constant or decaying SFR are not in¬ 
cluded in this table. ^ Mergers with burst intensities > 20. 
^ Intensity of the current SF burst with respect to the previ¬ 
ous average SFR. ^ Age of the current burst. ^ Fraction of 
the stellar mass formed in the current SF burst. 


they obtain a median f^ of 0.1, while ranges between 
0.6 and 0.9 in our sample, indicating that, in contrast to 
ULIRGs, dust with T < 50 K dominates the IR SEDs of 
LIRGs. 


The res olved stellar populations of lo cal LIRGs were ex¬ 
amined by Alonso-Herrero et al. (2010) for nine objects in 
our sample using optical integral field spectroscopy. These 
nine objects are mostly spirals or weakly interacting galax¬ 
ies. They found a higher contribution of old stellar pop¬ 
ulations to the total stellar mass in these LIRGs than in 
ULIRGs. This agrees with our result of lower mass frac¬ 
tions formed during the current burst of SF in spiral-like 
LIRGs than in merging LIRGs. 


5.2. Ha vs. IR luminosity. Tracing the last 100 Myr SFH 


Ha and IR luminosities are two fundamental SFR tracers. 
The former traces the most recent SFR since its luminosity 
is proportional to the number of ionizing photons, although 
it is necessary to apply an extinction correction that is not 
always accurately known. The latter, however, traces the 
obscured SFR, and it is virtually unaffected by foreground 
extinction, although it depends on the extinction level in 
the star-forming regi ons (e.g., objects with no extinction do 
not emit in the IR; |Ke nnicutt fc Evans| [2012). 


Figure [6] shows that our calibrations between the IR lu¬ 
minosity and the SFR at 1-10 Myr and 10-100 Myr, both 
^(9.64 ± 0.18 ) L 0 /(M 0 y r -1 ), are in very good agree¬ 
ment with the Kennicutt (jl998) calibration of 9.57 L 0 / 


(M 0 yr x ) (after correcting for the different IMF assumed). 


We also show in this figure the relation between the extinc¬ 
tion corrected Ha luminosities derived from the SED anal¬ 
ysis with the SFR. The logL corr (Ha)/SFR ratio is (7.14 
± 0.09) L 0 /(M 0 yr 1 ), so similar, although slightly lower 
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b - SFR i_ 10 Myr (M 0 yr ) b - SFR 10 _ 10 o Myr (M© Y r ) 




A v A v 

Fig. 6. IR and Ha luminosity/SFR ratio as a function of 
the fitted total extinction, A v considering the 1-10 and 10- 
100 Myr to average the SFR. The median ratio (solid red 
line) and the scatter (dashed red line) are indicated in each 
plot. The contribution of the AGN torus has been sub¬ 
tracted when present. The color of the symbols is as in 
Figure [4j 


than the calibration given by Kennicutt (1998) for the H a 
luminosity, 7.33 L 0 / (M e yr _1 p 

In the same way, we obtained the SFR calibration for 
the uncorrected Ha luminosity. The correlation between 
the L 0 bs(H<a)/SFR ratio and the extinction (A v ) is clear. 
However, since the average integrated extinction of these 
LIRGs is n ot extremely high, r v =0.5-3 or A v =0.6-3.3 mag 


(see also (Alonso-Herrero et al., 2006|)), the observed Ha 
luminosity could be used to estimate the integrated SFR in 
similar galaxies within a factor of ^3. 

We find that both the Lir/SFR and the L corr (Ha)/SFR 
ratios are approximately constant and do not clearly de¬ 
pend on the age interval used to average the SFR. However, 
part of the scatter of these relations can be explained by 
the differences in the SFH and the extinction level. This 
is clearer in Figure [7] where we plot the Air/L corr (Ha) ra¬ 
tio as a function of the extinction comparing the observed 
galaxies and our models predictions. We expect this ratio 
to depend mainly on the relative intensities of the current 
SF and the average SFR during the last 100 Myr because 
Ha emission is only produced by young stars, whereas older 
stars contribute to the IR luminosity too. But this ratio also 
depends on the extinction because the fraction of absorbed 
photons later re-emitted in the IR increases with the ex¬ 
tinction. The latter effect is apparent in this figure since 
the Tir/L corr (Ha) ratio increases, with A v tending asymp¬ 
totically to the ratio expected when 100% of the stellar light 
is absorbed and re-emitted by dust. 


In our sample, the age effect is more subtle because, 
according to the derived SFH, the SFR during the las t 
100 Myr has been approximately constant (see Section 5.1). 
In any case, to test this dependence we computed the ex¬ 
pected Tir/L corr (Ha) ratio as a function of the extinc¬ 
tion for several log SFRi_i 0 Myr/SFRio-iooMyr ratios. The 
model ratios are plotted in Figure [7] and we can see that our 
sample of LIRGs is, as expected, in agreement with the ra¬ 
tios for nearly constant SFH in the last 100 Myr. However, 
it should be noted that this ratio shows a strong dependence 
on the SFH. 

For comparison, the empirical relation between the cor¬ 
rected Ha luminosity and a linear combination of the ob- 
served Ha luminosity a nd the IR luminosity derived by 


Kennicutt et al. (2009) for a sample of nearby galaxies 
is plotted in Figure (7p The agreement between the data, 
our mo dels with log SFRi_ioM vr /SFRio_ioo Mv r = 0 — 1, 
and the |Kennicutt et al.| (|2009|) relation is good within the 
uncertainties. This comparison suggests that the empirical 
relation is only valid for galaxies with constant or nearly 
constant SFH during the last 100 Myr. Likewise, this dia¬ 
gram can be readily used to distinguish between galaxies 
with approximately constant SFH, such as those in the M-S, 
and galaxies with a rece nt (less than 100 Myr old ) SF burst 
that would lie below the Kennicutt et al. ( |2009 ) relation. 



A v 


Fig. 7. LiR/corrected Ha luminosity ratio as a function 
of the total extinction, A v . The dashed lines represent 
the expected ratio derived from our models for various 
logSFRi_ioMyr/SFRio-iooMyr values (red=2, orange=l, 
turquoise=0, blue=-l, and purple=-2). The dot-dashed 
black line is the ratio calculated from empirical relation 
WHa} = Lobs (Ha) + 0.0024 Lir derived by [Kennicutt| 


et ah I (2009) for a sample of nearby galaxies and the dotted 

ir 


black lines are the uncertainty of that calibration. The con¬ 
tribution of the AGN torus to the IR luminosity has been 
subtracted when present. The color of the symbols is as in 
Figure |4j 


6. Summary and conclusions 

We modeled the integrated SED of a sample of 38 local IR 
bright sources with IR luminosities ranging from log Lir/ 
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Lq = 10.2 to 11.8 and a median of logl/iR/L© = 11.0 
belonging to 29 systems with logLiR/L© = 11.0 — 11.8. 
The SEDs include new optical g, r, and H<a narrow-band 
imaging obtained with the NOT telescope combined with 
archival observations from UV to far-IR. To fit the SED, 


we modified the Bayesian code magphys (da Cunha et al. 
2008) and created a set of stellar population synthesis mod¬ 


els and dust models optimized for these objects. This SED 
fitting approach takes into account the balance between the 
energy absorbed in the UV/optical spectral ranges and that 
re-emitted in the IR. Except for three LIRGs that might 
host a deeply obscured energy source (AGN or SF), the 
SED models are able to reproduce the observed data well. 
The main results are summarized in the following: 

1. We classified the galaxies of our sample in three groups 
according to their SFH: objects with a recent burst 
of SF (47%), objects with a constant SFH (47%), and 
objects with decaying SFR (6%). In all cases the aver¬ 
aged SFR during the last 100 Myr seems to be relatively 
constant within a factor of three. 

2 . The intensity of the recent SF burst with respect to the 
previous averaged SFR varies between a factor of 2 and 
20 and it began ^30-300 Myr ago. The most extreme 
bursts (intensities > 10) are associated with mergers 
and high sSFR (> 0.8 Gyr -1 ). Similar to local ULIRGs, 
these objects would be compatible with a large part (up 
to 30%) of their current stellar mass being formed dur¬ 
ing the ongoing SF burst event. The rest of the stabursts 
include galaxies with varied morphologies: mergers, dis¬ 
turbed spirals, interacting pairs, and isolated disks. 

3. Galaxies with constant SFH in our sample tend to be 
more massive than those with a burst of SF (median 
log M*./ Mq = 11.0 vs. 10.6). Most of them lie just 
slightly above (<2 a) the M-S of galaxies, but they reach 
the LIRG luminosity threshold due to their high stellar 
masses. 

The calculated SFHs of the galaxies in our sample reveal 
that all of them had SFR and stellar masses in very 
good agreement with local M-S ~1 Gyr ago (z ~ 0.075). 
Their stellar masses were between 10 9/r and 10 11,4 M© 
and their SFR between 0.5 and 20M©yr _1 . It is likely 
that only one or two of them would have been classified 
as LIRGs if they would had observed 1 Gyr ago. 

5. We find that the L(IR)/L corr (Ha) vs. integrated op¬ 
tical extinction (A v ) relation derived for our galaxies 
is in goo d a greem ent with the empirical correlation of 


4. 


Kennicntt et al.| ( 2009| ). Our models show that this re¬ 
lation holds only if the SFR has been approximately 
constant during the last 100 Myr. Deviations from this 
relation can be used to identify galaxies with rapidly 
changing (increasing or decreasing) SFR during the last 
100 Myr period. 


Similar studies covering wider IR luminosity ranges and 
including integrated spectroscopic information will be cru¬ 
cial for obtaining a more detailed evolutionary history of 
local U/LIRGs, as well as for studying the differences and 
similarities with their high -z counterparts. 
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Appendix A: Optical spectra of nine galaxies 

We present the optical spectra of nine galaxies in our 
sample without a previously published [Nil]/Ha ratio, to 
best of our knowledge. Their reduced spectra were avail¬ 
able through the Smithsonian Astrophysical Observatory 
Telescope Data Center. They were obtained between 1998 


and 2006 with the FAST Spectrograph (Tokarz & Roll 
1997) on the Mount Hopkins Tillinghast 1.5 m reflector. 
The slit width was 3" and the spectra cover the spectral 
range between 3700 and 7500 A with a dispersion of 1.5 A 
per pixel. The integration times were between 600 and 
1500 s. The spectra are not flux calibrated but they can 
be used to measure line ratios between transitions close in 
wavelength. 

In the spectra shown in Figur jAT] we measured the 
fluxes of the H/3, [O ill] A5007A, Ha, and [N il] A6584A emis¬ 
sion lines using a single component Gaussian fit, except 
for CGCG 468-002 NED01, which shows a blue-shifted 
broad Ha component. The FWHM of the nar row co mpo¬ 
nent in this galaxy is ^500 km s -1 (Alonso-Herrero et al. , 
2013), whereas the broad Ha component has a FWHM 
of 3100T190 kms -1 and is blue-shifted by 180±60kms -1 . 
The observed line ratios are listed in Table O H[3 is not 
corrected for stellar absorption. 

We used the standard op tical d i agnos tic diagram [N n] / 
Ha vs. [Om]/H/? (Baldwin et al. 1981) to determine the 
nuclear activity classification. We used the boundary lim¬ 
its between Hu, composite an d AG N galaxies proposed by 
Kewley et al. (2006). Figure A.2 shows that four of the 
galaxies lie in tne composite region of the diagram, one in 
the Hu region, although close to the Hil-composite bor¬ 
der, and one galaxy, CGCG 468-002 NED01, is classified 
as AGN. Since in this object we only detect a broad com¬ 
ponent in Ha we classify it as Syl.9. For UGC 03405 and 
NGC 7753, the H/? and [Oiii]A 5007A transitions were not 
detected so we excluded these objects from the diagram. 
However, the high [N ll]/Ha ratio in these two sources, to¬ 
gether with the absence of [O III] detections, which is bright 
in AGNs, suggests that these are composite galaxies. 


Table A.l. Optical line ratios and nuclear activity classi¬ 
fication 


Name 

[Omj/H/3 

[Nil]/Ha 

Class. 

NGC 876 

0.75 

0.57 

composite 

NGC 2389 

0.45 

0.53 

composite 

CGCG 468-002 NED01* 

21.7 

0.98 

Syl.9 

CGCG 468-002 NED02 

0.89 

0.51 

composite 

UGC 03405 


0.75 

composite? 

UGC 03410 

0.46 

0.45 

H II 

NGC 7753 


1.02 

composite? 

NGC 7752 

0.60 

0.33 

composite 

NGC 7769 

0.91 

0.56 

composite 


Notes. Line ratios corresponding to the narrow component. 


Appendix B: Likelihood with detections and upper 
limits 

In this appendix we briefly describe how the up per lim - 
its are included in our Bayesian analysis (see, e.g, Gregory 
2005| |Bohm fc Zech||20l0| for a general description of the 



Fig. A. 2 . [Nn]A6584A/Ha vs. [O iii]A 5007A/H^ diagnos¬ 
tic diagram for the nuclear spectra of six LIRGs. The solid 
and dashed black lines mark the empirical sepa ratio n be¬ 
tween Hu, composite, and AGN galaxies from|Kewley et al. 
(2006). 


Bayesian approach). We let Fi and Oi be the flux and la un¬ 
certainty measured for a galaxy in the band i. If the galaxy 
is not detected, we measure the nai upper limit. Likewise, 
the prediction of the model k for the flux of band i is Mi(k). 
The likelihood is defined as 

L k = ~[[f i (F i ,a i ,M i (k )) (B.l) 


For the detections we assume that fi follows a normal 
distribution 


f?(F i ,a i ,M i ) = 


VTt, 


: exp 


(Fi - Mjf 

2a? 


(B.2) 


On the other hand, to obtain the likelihood for the upper 
limits, / z n , we first calculate the probability for the flux to 
have an arbitrary value Ri. The unknown background in 
the aperture used to measure the flux is Bi. The standard 
deviation of the background is cr*, and for simplicity we 
assume that the mean background of the image is zero. If 
the galaxy flux is Ri but we do not detect it on the image 
at a nai level is because Ri + Bi < nai. If the background 
follows a normal distribution the probability of this is 

rnn a, Ri) = $ ( B -3) 

where <f> is the cumulative distribution function of the stan¬ 
dard normal distribution. Therefore the likelihood value for 
the upper limits is 


f?{F i ,<r i ,R i ) = 


1 + erf 


nai ~ Ri 
<?i \[2 


(B.4) 


Then whe n substituting Equations |B.3| and |B.4| into 
Equation |B.1| 
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480 485 490 495 500 505 640 650 660 

rest wavelength (nm) 



Fig. A.l. Nuclear optical spectra of nine galaxies of our sample obtained with the FAST Spectrograph in the ranges 
480-510A and 637-680A. The dashed lines mark the position of the H/?, [O ill] A5007A, H<a, and [N ii]A 6584A transitions. 


£ *=n 

H 

n 


\f2nfai 


: exp 


(F h - M h (k)) 2 
2a ? 


1 + erf 


n<jj 2 - Mj 2 (fc) 
&i 2 V2 


(B.5) 


where i\ and i 2 are the subindices for the detections and 
non-detections respectively. The logarithm of Equation |B.5| 
is 


in L k =yy 

h 

E ln 


(F h - M n (k)y 


2(7? 


1 + erf 


+ 


(B.6) 


na i2 - M i2 (k) 


3 V2 


+ C 


= - ^xI + E 111 


1 + erf 


nail - M i2 (k) 
&i 2 V2 


C 


Thus, we assign the following probability to model k in 
the parameter inference process 


L k = exp f-Rfe) hi 


1 + erf 


ncr i2 - Mj 2 (k) 
Vi 2 \f2 


(B.7) 
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Appendix C: Multi-wavelength imaging of the LIRGs 


MCG +12-02-001 



NOT g ' ■ 


50 25 0 -25 -50 

RA (arcsec) 

NGC 876 


-50 ■ 


MIPS 24pm 


PACS 70pm 


PACS 160pm 


... ;i„ . 

: NOT Ha+[NII] -r 2MASS K 1 IRAC 3.6am IRAC 8.0|jm 


- SPIRE 250gm 





2MASS K 


IRAC 3.6pm 


MIPS 24pm 


PACS 70pm 


PACS 160pm 


-25 ■ 


IRAC 8.0pm 


SPIRE 250pm 


50 25 0 -25 -50 

RA (arcsec) 


Fig. c.i. Same as Figure [l] 
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Fig. c.i. Continued. 
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Fig. c.i. Continued. 
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Fig. c.i. Continued. 
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Fig. C.l. Continued. 
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Fig. c.i. Continued. 
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Fig. c.i. Continued. 
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Fig. C.l. Continued. 
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Fig. c.i. Continued. 
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Fig. c.i. Continued. 
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Fig. C.l. Continued. 
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Fig. c.i. Continued. 
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Fig. C.l. Continued. 
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Appendix D: Best fitting model results to the SED 
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